Phase diagram of hole doped two-leg Cu-0 ladders 
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In the weak coupling limit, we establish the phase diagram of a two-leg ladder with a unit cell 
containing both Cu and O atoms, as a function of doping. We use bosonization and design a specific 
RG procedure to handle the additional degrees of freedom. Significant differences are found with the 
single orbital case; for purely repulsive interactions, a completely massless quantum critical region 
is obtained at intermediate carrier concentrations (well inside the bands) where the ground state 
consists of an incommensurate pattern of orbital currents plus a spin density wave (SDW) structure. 



The challenging physics of strongly correlated systems 
provides a unique opportunity to test many proposals for 
new, unconventional quantum states of matter. In that 
respect, ladders constitute a particularly interesting case 
[H . These ID systems behave quite differently from single 
chains. One can show for instance that, for purely repul- 
sive interactions, they favor superconductivity in their 
ground state. Understanding their properties - both ex- 
perimentally 0] and theoretically - is thus interesting in 
its own right but also could help us gain valuable insight 
into the elusive physics of the two dimensional cuprate 
superconductors. 

Most studies of ladder compounds model these systems 
with a single orbital per unit cell [1, !3| . Using a renormal- 
ization group (RG) analysis of the Hamiltonian expressed 
in bosonic variables, a phase diagram can be derived in 
the weak coupling limit. As pointed out [3=|, the param- 
eter that may be safely tuned to arbitrary values in the 
weak-coupling limit is the doping d and its variation pro- 
duces a sequence of states, labelled CnSm, with n (m) 
gapless charge (spin) modes. For repulsive on-site Hub- 
bard U terms, a CI 50 d-wave "superconducting" phase 
is found away from half filling {5 ^ 0). Relaxing the 
constraint on the magnitude and on the sign of the inter- 
actions and extending their range to more distant sites 
allows one to promote other types of orders such as or- 
bital antiferromagnetism (OAF) [5j. This state is also 
known as a flux phase and was examined in the context 
of the two dimensional Hubbard model [1, 0| • 

The question of whether orbital currents could exist in 
cuprate materials has received much attention. Analyt- 
ical ^ and numerical studies of a three band model 
of the copper-oxygen planes predict that, in the large U 
limit, a strong Coulomb repulsion Vcu-o between near- 
est neighbor Cu and O atoms favors such a phase. Recent 
experimental data seem to support that picture [lo| , but 
more studies are needed to confirm this scenario. 

We thus revisit models for ladders and include the oxy- 
gen atoms. Here, we focus on the issue of whether their 
presence causes any significant changes [lH and in par- 
ticular whether orbital phases might exist for realistic 
choices of microscopic interactions, ior S ^ 1^. We 
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FIG. 1: (Color online) (a) Molecular structure of the unit cell 
showing the hopping and interaction parameters included in 
the hamiltonian (b) Orbital current pattern and SDW modu- 
lation (bold arrows) in the C2S2 phase. The chain direction 
is vertical, kpo is the Fermi wavevector in the band. 



establish, in the weak coupling regime, the phase di- 
agram as a function of hole doping of two-leg ladders 
whose unit cell contains both Cu and O atoms, with 
on-site repulsions Ucu (Uo) on the Cu (O) sites and 
a nearest-neighbor Vcu-o Coulomb term. We use RG 
techniques to map out the flows for the bosonized ver- 
sion of the model. In contrast with the case of a single 
orbital ladder we find that for an intermediate range of 
dopings 6ci < S < Sc2 a fully massless phase is stabilized. 
The value of (5ci(c2) depends on the bare magnitude of 
the Hubbard terms and/or the interoxygen hopping tpp 
which we treat as tunable parameters. Furthermore, in- 
creasing tpp beyond a minimum value t™" promotes, for 
all 5 > Sci, an incommensurate orbital current state. It 
is similar to one of the patterns advocated by Varma 
(see Fig. [TJ)) . The corresponding phase has an additional 
SDW (CDW) character for S < Sc2 {S > Sc2)- 

We consider the two leg ladder of Fig. [Iji where the 
relevant hopping parameters and interactions are shown. 
Ht is the sum of contributions describing carrier hops 
plus a term proportional to e = Ecu — Eq, the differ- 
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ence between the Cu and O site energies. Relevant val- 
ues pertaining to selected copper oxide ladders have been 
computed in LDA and we use here = 1, e = 0.5 
in units of t. Hx is diagonalized in momentum space, 
and since e and the various t's are of the same order, one 
can safely neglect the high energy orbitals. We are left 
with two low lying bands crossing the Fermi energy Ep, 
one bonding (0) and one antibonding (tt) combination of 
chain states. The energy dispersion is linearized near Ep . 
From the fermionic densities, we introduce charge (c) 
and spin boson (s) fields (j) for each specie. Ht is diago- 
nal when expressed in terms of operators /i = c or s, 
J' = or TT in the Bqtt basis. The non linear terms of Hint , 
denoted by Hj^^, have a simple form in the basis, 



where 



int 

1/V2(0, 



When 



density-density interactions are included, the quadratic 
part of H is diagonal in the Bo basis where we introduce 
operators with A = 1,2,3,4 (1,2 are s modes, 3,4 are 



c modes): 

^0 = E/ ^[{^xK.Knn.r + i^Jid^M'] (1) 

The matrix S which defines the rotation of the Bo basis 
with respect to is given by 



f Pi Qi \ 

J_ -Qi Fi 

V2 P2 Q2 

V -Q2 P2 J 



(2) 



Pi and Qi are expressed in terms of angles a for the spin 
part and (3 for the charge part; Pi{2) = cosQ:(/3)-|-sina(/3) 

and Qi(2) — cos a{/3) — sma{(3). In the B^ basis H 

reads: 



NL 
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H. 
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gic / dr cos(20s+) • cos(26'c-) + 5ia / cos(20s+) • cos(20s-) — g2c / dr cos{29c~) ■ cos{2(l>s-)+ (3) 



54a / dr cos(2(/)s_)-cos(20s_)+5i j (ircos(2(/)s+)-cos(2(/)s_)+g2 / dr 8iTi{2(j)s-) sm(2(l)s^)+g\\c j dr cos{29c-)-cos{26s^) 

r 



Here, we use the same convention for the Klein factors 
as in [l5l |. Subscripts 1 to 4 have the standard g-ology 
meaning and labels a to d refer to processes involving the 
and/or tt bands. The two gid terms, e.g, describe events 
where one right- and one left- moving fermion, both be- 
longing to the same ( or tt) band, backscatter within 
that band, gi and 52 correspond to the sum and to the 
difference of these "ld"-type processes respectively, and 
52 7^ when the O atoms are included. The g^a term has 
a non-zero conformal spin so that additional interactions 
Gp(t) ~ cos(40s_ (6's_)) are generated during the flow. 
They are included in our calculations. Since we are con- 
cerned with a priori incommensurate values of S we drop 
all umklapp terms. We renormalize the couplings in ([3|) 
following the usual RG procedure, where one integrates 
out high energy states. This sequence is straightforward 
when the quadratic part ([T]) is expressed in the Bo basis, 
since one deals with simple gaussian integrals but when 
we express in the Bo basis this involves Pi and Qi 

coefficients. Each RG step then generates cross-terms in 
Hq, which implies a rotation of Bq with respect to the 

_B_i basis. It is thus important to include the change in 

S during the flow. After the n*'* iteration, we denote by 

(a„,/3„) the angles between Bo and B^ and by K^'^\ 

the parameters in the Bo basis. We perform the (n-l- 1)*'' 
RG step in the Bo basis, which changes K^"'\ (see (P)) 



and introduces cross-terms. We apply 5'~^(a„,/?„) to 
ffo, which takes us back to the (fixed) B^ basis. Fi- 
nally we determine the new angles (a„+i,/3„+i) which 
are required to make Ho diagonal again, with new pa- 
rameters i4r("+-^^A- 

Proceeding in incremental steps gives the additional 
RG equations for the rotation of the Bo basis 



dcot2a{l3) _ {{dK 1(^3,) - dK 2{i))t&nAa + dB i2(zi)) 



dl -^1(3) - K2(A) 

The equations for the off-diagonal terms c?-Bi2(34) are 



■dl-^ 
(4) 



dBri 
dl 



PiQiiiala 



Gl)-K^K2{9la+glc 



+gl + GD) - KiK2hiPi)gig2 
dBai ^ , 2 



dl 



P2Q2{gt 



I) 



(5) 



with h{Pi) = ((FiQi)2 + 0.25(^1^ _ Ql))-\ Details of 
the RG equations for the various g, K and for the ratio 
of the Fermi velocities in the and tt bands will be given 
in a forthcoming publication [I61] . 

Using the above equations we establish the phase dia- 
gram for the ladder. In agreement with ^ we find that 
a = ^''2vJ''^ ratio of Fermi velocities in the and 
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TT bands) controls the behavior of the differential equa- 
tion system. When -j- is constant, a depends only on 
d. Two main factors may significantly affect the phase 
diagram that was predicted for two leg Hubbard ladders 
with a single orbital per site: one is the asymmetry in 
the g terms due to the fact that the projections of the 
Cu and O orbitals onto the and tt bands have unequal 
amplitudes and one is the influence of the extra parame- 
ters Uo, Vcu~o and tpp. We first investigate the impact 
of the asymmetry, so we set Uo — Vcu-o — tpp — 
and we choose small initial values for Ucu (in the range 
10~6 _ 10-1). 

(a) Small doping range. For small 6 (a), cot 2a 
and cot 2/3 —> as the flow converges towards the fixed 
point, thus Bo In this case, g2, g4a, Gp, Gt are ir- 
relevant while 9c- and are ordered (0c_ — 0, (j)s+ — 
mod 27r). This is the CISO phase [3| where only the 
total (-I-) charge mode is massless. For the s— (spin- 
transverse) mode, terms involving the canonically con- 
jugated fields 4>s- and 9s- are relevant and competing, 
d-type superconducting fluctuations (SCd) dominate if 
(j)s- is locked at 0, while OAF is favored if 9s- = 0. 
Here, SCd is always more stable for repulsive Ucu- This 
property holds only for 5 < Sd — 0.2, where the spin and 
mass gaps go to zero. 

(b) Large doping range. For 6 > dc2 = 0.28, cot 2a 
and cot 2/3 oo (with opposite signs), so Bo Sott- 
In this regime only gi ~ —g2 are relevant and they lead to 
a state with one massive spin mode (in the band) . This 
is the C251 phase. The slowest decay of correlations is 
observed for the CDW operator in the band. Fluctua- 
tions in the tt band favor a SDW state (when logarithmic 
corrections due to a marginal operator proportional to 
51+52 > are included) but they are subdominant. If 
6 is just above Sc2, gi and 172 increase very slowly during 
the flow and one needs to choose larger values for the 
bare Ucu (still much smaller than t) to reach the asymp- 
totic regime with a gap in the spin mode. In contrast 
with the case of a single orbital per site, the C2S1 phase 
is stable, even for dopings such that Ep is close to the 
bottom of the tt band where a is very large (in that limit, 
we cannot linearize the energy spectrum, but we use di- 
agrammatic techniques [^). This comes from the fact 
that for unit cells with Cu atoms only, g2 is accidentally 
equal to zero. When O atoms are included (or when 
Vcu-o 7^ 0) the initial 92 is non-zero and 52 is always 
relevant. The nature of this C2S1 phase is discussed in 
the next paragraph. 

(c) Massless regime in the (Sd, 602) range. As S ap- 
proaches the critical end points Sci and 5c2 respectively 
from below (in the CISO phase) and from above (in the 
C2S1 phase), gaps in the spin and/or in the charge sec- 
tors go to zero. All spin and charge modes are mass- 
less in the entire range of dopings S G (0.2; 0.28). ^ 
and ^ are very large and the fixed point values of /3 
(a) just below and just above 6ci (<Jc2) are significantly 



different. So we approach the critical points from the 
massive phases at both ends; we discard couplings which 
flow to zero and thus obtain a simpler set of equations. 
Next we single out terms in (|3][S|) which produce large 
values of the derivatives in this range and determine the 
fixed point value of cot 2a(/3). It allows us to write down 
a minimal set of differential equations for the couplings 
and to determine those which are relevant in the doping 
range {6ci, Sc2). We first consider dopings close to Sc2- 
This corresponds to an initial value of cot 2a equal to 
one. The signs of {dKi — dK2) and dBi2 are the same 
and positive whereas the sign of [Ki — K2) is negative 
so that, according to ([3]), |cot2a| decreases to zero be- 
low 5c2 while above it increases to infinity. Below 5^2 1 
gi and 52 are not relevant and the system flows to the 
C2S2 phase while above they are relevant, leading to the 
C2S1 phase. For 5 close to 6ci = 0.2, (di^3 - dKA) and 
(ii?34 have opposite signs. Depending on which of the 
two terms dominates, | cot 2(3\ goes to zero or to infinity. 
At 5ci they are exactly equal. | cot 2/3| 00, for 5 > 5ci, 
but, since 0.5 < Ka < 1, one finds that all interband 
couplings (as well as higher order cos(60c) terms with 
b — 4,6, ..) are irrelevant. For 5 G (0.2; 0.28), all interac- 
tion terms are irrelevant and Bott {B-\ ) is the eigenbasis 

for the charge (spin) modes. A numerical solution of the 
full set of RG equations confirms this statement. The 
existence of this massless regime is essential to maintain 
spin rotational symmetry in this doping range. The na- 
ture of the C2S2 phase can be determined in the frame- 
work of the Luttinger liquid description. In that phase, 
K4 (K^) which corresponds to the (tt) band is smaller 
than (around) one. Dominant fluctuations occur in the 
band, and this case is equivalent to treating a single 
chain problem. The only marginal couplings are gi and 
92 (51 > 52)- Including logarithmic corrections allows 
us to identify the slowest decaying correlation function 
and we find that in the {Sd, Sc2) doping range, a SDW 
in the band {SDW{o)) is favored. In the C2S1 state, 
gi ~ —g2 < are relevant which gives a gap in the spin 
sector of the mode. In that regime, fluctuations in the 
band dominate, and one finds that the CDW{o) state 
is the slowest decaying one. 

Next, we "turn on" the parameters Uq, Vcu-o and tpp 
and we examine their impact on the phase diagram. In 
the doping range covered by case (a), SCd becomes less 
dominant over OAF when we increase the (positive) bare 
value of Uo or Vcu-o at fixed tpp but it is still the phase 
with the lowest free energy. One would need to assume 
a very large attractive bare Vcu-o to cause a transition 
12 1 to an s-type SC phase (^s- — 0, 4>s+ — 0, 9c- = 7r/2) 
15| which persists even as Ep approaches the bottom of 



the TT band. As far as case (b) is concerned, we observe 
a reduction in the size of the gap for positive Uq and 
Vcu-o, while for very large attractive Vcu-o the s-SC 
phase re-enters. In case (c), increasing Vcu-o has little 
effect on 5c2 but it shifts 5ci to higher values. An un- 
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FIG. 2: (Color online) Knight shift on the outer O sites calcu- 
lated for different phases: solid line CISO, dashed line C2S2, 
dash dotted line C2S1. 



physically large ratio Vcu-o/Ucu ~ 5 would be required 
to suppress the massless phase and to observe a reen- 
trant CISO phase with superconducting fluctuations so 
that in the relevant case Vcu~o < Ucu the intermediate 
massless phase does exist. 

The interoxygen hopping has a more significant effect. 
Increasing the value of tpp causes a concomitant decrease 
of 5ci and Sc2- For tpp = 0.5 -a value pertaining to 
Cu — O ladders Il3|- their values are about half that 



quoted for tpp = 0. If tpp > i™", a new phase is favored 
when S > 6ci- This state has both orbital current and 
DW fluctuations {DW = SDW (CDW) for the C2S2 
{C2S1) regime) and it shows similarities with one of the 
patterns advocated by Varma Q (see Fig. [1)3). This 
current phase is an eigenstate of the Bott basis (it is in- 
variant under the exchange of the two legs) in contrast 
with the other Varma pattern or with the usual OAF. 
The pattern has an incommensurate spatial periodicity 
~ kpQ. The amplitude of this order parameter is a sum 
of current operators between hnks of the Cu-O loops, of 
the form tijIm{X*QXjo), where tij is the hopping param- 
eter from site i to site j within the same unit cell and 
Aio is the overlap between the {Cu or O) wavefunction 
at site i and the band eigenfunction. These quantities 
are of order one [l6j and change by only a few percent 
when S increases from 5ci to the bottom of the band. Due 
to current conservation, the weakest link between atoms 
determines the maximal current, and we find that, for 
i™" « 0.3, the "Varma" state dominates the DW(o). 
The presence of the O sites insures that Im{X*^\ja) 7^ 
(a = 0,7r); Otherwise, the current operator between Cu 
atoms has the usual interband form: cJ^c,ro-- 

Our predictions could be tested by performing NMR 
measurements on the Sri^-xC axCu24:04:i compound [l| 
where the hole content can be somewhat varied, as we 
find very different responses of the the spin modes for 



doping we find an activated behavior of K (and also T^^) 
and K{T = 0) = 0; for high doping the temperature 
dependance is similar but i^(r = 0) 7^ 0; finally, at in- 
termediate dopings the usual high temperature behavior 
if ~ T° (saturation to the LL value) extends to T — > 0. 
For instance, Fig. [5] shows the Knight shift predicted for 
the "outer" O sites (i.e interladder bridges) 

In conclusion, we have developed a new RG method to 
handle correlation effects in the weak coupling limit for 
two leg Hubbard ladders at generic filling, when oxygen 
atoms are included in the unit cell. We have found a 
ground state phase diagram where the CISO and C2S1 
phases are present at small and large dopings, as for the 
single orbital ladder, but also a new intermediate phase 
C2S2 which is completely massless. Dominant fluctua- 
tions in the C2S2 and C2S1 states correspond to orbital 
currents preserving the mirror symmetry of the ladder 
structure on top of a SDW(o) and CDW(o) respectively. 
The stability of this new phase to an interladder coupling 
and/or to large values of the bare U are under current 
investigation. 
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the Cu and O sites [16|. Knight shifts K and relaxation 



rates ^ can be evaluated in all three regimes. For low 
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